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Abstract. We obtained Hq diffraction limited data of 
the LBV star P Cyg using the ONERA Adaptive Optics 
(AO) facility BOA at the OHP 1.52m telescope on Octo- 
ber 1997. Taking P Cyg and the reference star 59 Cyg AO 
long exposures we find that P Cyg clearly exhibits a large 
and diffuse intensity distribution compared to the 59 Cyg's 
point-like source. A deconvolution of P Cyg using 59 Cyg 
as the Point Spread Function was performed by means 
of the Richardson-Lucy algorithm. P Cyg clearly appears 
as an unresolved star surrounded by a clumped envelope. 
The reconstructed image of P Cyg is compared to simi- 
lar spatial resolution maps obtained from radio aperture 
synthesis imaging. We put independent constraints on the 
physics of P Cyg which agree well with radio results. We 
discuss future possibilities to constrain the wind structure 
of P Cyg by using multi-resolution imaging, coronagraphy 
and long baseline interferometry to trace back its evolu- 
tionary status. 

Key words: star: P Cyg, techniques: adaptive optics, de- 
convolution 



1. Introduction 

Among the galactic Luminous Blue Variables (LBV), the 
supergiant P Cyg (HD 193237, Blla-|-) is both a historical 
prototype since its famous XVIIth century eruptions, and 
an intriguing paradigm in many aspects. In their exhaus- 
tive study of P Cyg's fundamental parameters, mas s loss 
physics and evolution, Nota et al. ( Nota et al. 1995 ) have 
held this star being the unique case of the so-called " Pecu- 
liar Nebulae" class LBVs. Indeed the 22 arcseconds (arc- 
sec) faint and spherical nebulae around P Cyg were only 
recently discovered by Barlow et al. ( Barlow et al. 1994 ) 
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who estimate its mass as 0.01 solar mass. This value sig- 
nificantly differs from those of more regular LBVs, ranging 
from IMq to A.2Mq (e.g. for AG Car) which also present 
more asymmetric nebula. 
Evolutive tracks suggest an initial mass of 48±6 Mq, 
and it's present mass is estimated to be at the most 
40 Mq, but lower masses (~ 30 Mq) are also reported 
in literature ([Lamers et al. 1983 , Lamers et al. 1985, and 
Turner et al. 199S| ). However, the fine spatial structure 
of this large amount of excreted matter remains to 
be detailed. P Cyg's relative proximity (^ 1.8 kpc. 



Lamers et al. 1983) represents an opportunity to observe 



its radiatively driven mass loss from the starting point 
out to the interstellar medium. Indeed, at 1.8 kpc, the 
central star radius (76it!0) corresponds to a tiny angle of 
0.2 milliarcsecond (mas) but the Hq emitting region ex- 
tends over several tens of arcsec, and radio emission seems 
to attain at even larger scales (Meaburn et al. 1999). 



The optical and radio observations depict so far an es- 
s entially clumpy d i stribution of ma tter both at large 
(Taylor et al. 1991, Nota et al. 1995) and small scales 



( Skinner et al. 1998 , Vakili et al. 1997 ), with temporally 
variable emission ( Skinner et al. 1997 ). These imaging ob- 
servations remain sparse and can loosely constrain the spa- 
tial and/or temporal evolution of the clumps in the nebula. 
Moreover, in the optical wavelengths, the star to the en- 
velope brightness ratio remains an obstacle for studying 
the central star immediate environment. In this paper we 
report an attempt to observe the Hq circumstellar envi- 
ronment of P Cyg during an AO run at OHP observatory 
on October 1997 using short exposures collected with a 
photon-counting camera. 

The paper is organized as follows. In section 2, we describe 
the AO-|-imaging instrumentation used for this study as 
well as data reduction and calibration procedures. Section 
3 describes the image reconstruction in which particular 
efforts have been made to test the validity of the PSF. 
The last section makes a critical discussion of our found- 
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Transrnited spectrL 



Fig. 1. Hq long exposures of P Cyg (top) and the refer- 
ence star 59 Cyg (bottom) recorded with the CP20 cam- 
era at OHP during October 4th 1997 using the AO sys- 
tem BOA from ONERA through a filter (AA=10nm). The 
images are smoothed using a moving average window of 
4x4 pixels. Note that the theoretical angular resolution of 
the 1.52m telescope at 65631 (1.22A/D=0.110 arcsec) is 
closely reached (broken first Airy ring is visible on 59 Cyg) 
whilst for P Cyg a more complex and resolved structure 
can be suspected. 

ing, and attention is given on the perspectives opened by 
the emergencies of new high angular resolution techniques. 
A need arise for a coordinated campaign for a global ap- 
proach of P Cyg, and LBVs environment. 

2. Journal of observations 

2.1. Instrumentation 



P Cyg was observed, in the context of a dark-speckle run 
( iBoccaletti et al. 1998D . The AO bench of the ONERA 
(Office National d'Etudes et de Recherches Aerospatiales) 
was installed at the Coude focus of the 1.52m telescope. 
The 90 actuators of BOA and its 1 kHz closed loop band- 
width enable a compensation of atmospheric turbulence in 



visible light (Conan et al. 1998) and provide typical Strchl 
ratio of 10% to 30% depending on seeing conditions. 
The restored wavefront feeds a dark-speckle optics de- 
signed to detect faint companions around bright objects. 
The beam is then focused onto the detector with an //976 
aperture, giving a fine pixel sampling of 144 pixels /arcsec 
or O.OOlarcsec/pixel (this oversampling of the images was 



: I ; 



6600 
Wovclcnght (ongstrom) 



Fig. 2. Synthetic representation of the Ha region of P Cyg 
multiplied by the filter transmission used for the adaptive 
optics imaging of P Cyg. The same filter was used for 59 
Cyg. About 45.6% of the transmitted light originates from 
the line continuum. 



imposed for the dark-speckle experiment) . The detector is 
a cooled CP20 photon-counting camera ( Abe et al. 1998 ) 
allowing single photon detection with a very low dark 
count (< 5. 10^^ photon /pix/s). The quantum efficiency 
of CP20 is less than 10% at 700nm. 

Near diffraction-limited images of P Cyg (to„ = 4.81, 
S„ = Blla+) and 59 Cyg (m^ = 4.74, Sp = B1.5V) 



have been obtained (Fig. m using a broadband H^ filter 
(Ao = 6563A, AA = lOOA) (Fig. |). 



2.2. Data analysis 

The data consist in two sequences of 20 ms short ex- 
posures, for respectively P Cyg, and the reference star 
59 Cyg of similar magnitude and spectral type. In a 
preliminary step, the short exposures recorded by CP20 
were cleaned for photon-centroiding electronic artifacts 
and then co-added to generate an equivalent long expo- 
sure of 99s for P Cyg and 370s for 59 Cyg. Despite the 
Coude configuration of the 1.52m telescope, rotation of 
the field between the two sequences was found negligible. 
We checked for non-linear effects introduced by the 
s o-called photon -centroiding hole and camera saturation 
( Thiebaut 1994 ). It appears that the incoming flux is 
much below the saturation limit, the double-photon oc- 
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Table 1. Turbulence characteristics derived from the 
wavefront sensor data of the adaptive optics system. At- 
mospheric parameters are given for the visible band. 



Star name 


PCyg 


59 Cyg 


Time 


20:30 TU 


21:23 TU 


Average intensity 

in sub-apertures (photons) 


51.567 


42.830 


ro (cm), open loop 
calculated with Lq 


5.83 


5.92 


ro (cm), open loop calculated 
with Lo and Zernike coef. 


~5.4 


~5.4 


I/O (m) 


2.96 


2.57 


Zernike coef. (2,3), cr^ {rad^) 


70 


40 


Zernike coef. (4,44), a'' (rad^) 


10 


9 



currence is absent. The total number of detected photons 
for the two images are very similar: 481563 for P Cyg and 
416309 for 59 Cyg. A coherent peak and a broken ring, 
featuring a triple coma aberration, are clearly visible, 
but 59 Cyg's image is sharper, indicating that P Cyg's 
envelope is possibly resolved. This is quite visible in 
Fourier space, as shown in Fig. |[ The envelope of P Cyg 
clearly appears as the well resolved low frequency part in 
the visibility curve. High angular information is obtained 
up to about 80 % of the cut-off frequency. The visibility 
curve is indeed clearly dominated by noise beyond that 
limit. 
Since 59 Cyg has been observed 53 minutes later than 
P Cyg, a variation of atmospheric conditions and/or 
a poor correction could be responsible of a noticeable 
change in the Point Spread Function (PSF) shape. From 
open/close loop data recorded by the wavefront sensor 
of BOA, seeing conditions have been derived by one of 
us (F.C.) : intensity in sub-apertures. Fried parameter 
evolution, noise evolution, Zernike coefficient variance 
evolution, etc... Some of these relevant parameters are 
summarized in Table 1 for both stars. In particular 
the Fried parameter (rg) exhibits no temporal variation 
and remains stable around 5.4cm between the 2 data 
sequences. Also, high order correction (Zernike coef. 
(4,44)) have similar variances in both cases, indicating 
similar adaptive mirror corrections for P Cyg and 59 Cyg. 
The largest difference appears for the tip-tilt corrections 
(Zernike coef. (2,3) of the Table 1). However we confirm 
that during the observations no saturation was detected, 
and we expect these modes be correctly compensated for 
both stars. Furthermore, parameters such as turbulent 
layers altitude, mean isoplanetic angle and wind velocity 
indicate a very good stability of the atmospheric condi- 
tions during the two records. 

Moreover, the photon flux stability can be checked on 
the short exposures statistics. For 59 Cyg, the flux is 
stable with a perfect Poisson variance. P Cyg's behavior is 
roughly the same, except some localized empty exposures. 




Fig. 3. Representation in contour plots of the modulus of 
the Fourier transform of P Cyg (top), 59 Cyg (middle) 
and their ratio, the visibility curve of P Cyg (bottom). 
The zero frequency is at the center of the images. The 
representations are in logarithmic scale. The circle gives 
the theoretical cut-off frequency of the telescope. 



due to electronic saturation, which have been removed. 
This preliminary analysis permits to assume that 59 Cyg 
can be confidently considered as P Cyg's PSF and a 
successful deconvolution becomes possible. 



3. Image reconstruction 
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3.1. Deconvolution of the raw image of P Cyg 

A deconvolution of the image of P Cyg using 59 Cyg 
as the PSF has been made using basically the itera- 



tive Richardson-Lucy algorithm (RL) (Richardson 1972, 
Lucy 1974) to which some improvements have been ap- 
plied as we shall describe thereafter. To improve the re- 
sult, a pre-processing of the image was first performed. 
All images were smoothed using a moving average win- 
dow of 4x4 pixels. To get rid of the effects of the coron- 
agraphic mask (clearly visible in Fig. ^ at the right of P 
Cyg's image), only a central region of 0.9x0.9 arcsec of the 
smoothed images was conserved. To make the deconvolu- 
tion procedure faster, the number of points of the images 
was further reduced by a factor 16. The resulting images 
of 32x32 pixels remain correctly sampled (0.028 arcsec 
compared with the theoretical Shannon limit of 0.044 arc- 
sec for the 1.52m telescope operated at 656. 3nm). 
The RL algorithm was then applied to these pre-processed 
images using the well-known iterative procedure: 



x'=+i(r) =x'^(r)-/i(-r) 



yir) 



h{r) ® x^{r) 



(1) 



where ® is the convolution symbol, r denotes the two- 
dimensional spatial position, y[r) denotes the image of P 
Cyg, h{r) the image of 59 Cyg taken as the PSF, and 
where x^ {r) , x^'^^ {r) are the reconstructed images at the 
iteration k and fc 4- 1 . 

To limit the instability that appears in the solu- 
tion, due to the amplification of noise, we stop the 
iteration number by using a comparison between the 
Fourier Transform of the reconstructed object at the 
iteration k and the Fourier transform of the image 



reconstructed by a Wiener filter (Lanteri et al. 1998, 
Lanteri H., Soummer R. et al. 1999). The main difficulty 



is then shifted to a correct determination of the Wiener 
filter. The comparison lead to choose an iteration num- 
ber of 245. The result of the deconvolution is given both 
in gray levels and in contour plot in Fig. 0. To make the 
envelope clearly visible, the contour levels are not equally 
spaced. The general pattern is that of a bright star, not 
resolved by the 152 cm telescope, surrounded by an ex- 
tended envelope with bright spots. 

These same results were obtained when we used for the re- 
construction the RL algorithm regularized by a Tikhonov 
t erm, and more specially using the Laplacian operator 
(iLanteri H., Roche M. et al. 1999]) . The result of this de- 



convolution is presented in Fig. |5|for the iteration 250 with 
a regularization factor equal to 0.01. Neglecting the fine 
structure of the envelope, a normalized circular average of 
the image centered on the star may be fitted by the sum of 
two Gaussian curves, of the form: ae~''l''l -t- (1 — a)e~'' I*"' 
where a ~ 0.85, b ~ 618.2 and b' ~ 25.9 if we express 
r in arcsec. This corresponds to an equivalent width of 
the envelope of about 0.4 arcsec (-^). The radial curve is 
shown in a semi-logarithmic scale in Fig. |6[ In this model 




Fig. 4. Reconstructed image of P Cyg with the RL algo- 
rithm stopped at iteration 245 . Top: Representation in 
gray levels in a linear scale. Bottom: Representation in 
contour plot. Contours levels are not equally spaced and 
correspond to : 100, 80, 60, 40, 30, 15, 12.5, 9, 6.5, 5, 3.5, 
2, 1.2 percent of the image's maximum. North is at the top 
and East is at the right of the images. This corresponds 
to a rotation of 104.7° clockwise with regards to Fig. ^. 



the total integrated flux produced by the envelope (r = 
(l-a)b/ab') was found to be about 4 times larger than 
the one produced by the central star. This ratio is not an 
absolute parameter; it depends of course of the spectral 
bandwidth of the experiment (Fig. g) . 
We may conclude from this first analysis that the en- 
velope of P Cyg is well resolved by the central core of 
the PSF of the 1.52 m telescope, corrected by the adap- 
tive optics system. However, the envelope remains com- 
parable in size with the residual halo of the PSF. This 
assumes at least that we are free of anisoplanatism prob- 



lems (Fusco et al. 2000); this absence of variation allows 
the use of conventional deconvolution methods. 
The question then may arise whether some of the fine 
structures we discovered in this envelope may be due 
to artifacts of the reconstruction process. Among oth- 
ers problems, a very important one is a possible vari- 
ation of the PSF during the experiment (from P Cyg 
t o 59 Cyg), as discussed recently by Harder and Chelli 
( Harder and Chelli 2000 ). These authors show that a lo- 
cal non-stationary turbulence may produce strong residual 
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Fig. 5. Reconstructed image of P Cyg with the RL al- 
gorithm regularized by the Laplacian operator. Contours 
levels are not equally spaced and correspond to: 100, 80, 
60, 40, 30, 15, 12.5, 9, 6.5, 5, 3.5, 2, 1.2 percent of the 
image's maximum. North is at the top and East is at the 
right of the image. 




Fig. 6. Thick line: Radial profile of the object recon- 
structed by RL stopped at iteration 245 (image normalized 
to one). Dashed lines: Gaussian curves. Thin line: Sum of 
the two Gaussian curves. Intensity is in logarithmic scale. 



aberrations (clearly visible in the first diffraction ring of 
their figure 17). At worst, can we imagine that the ob- 
served structures in our image result only from variations 
in the PSF, from P Cyg to 59 Cyg? 

A comparison between the reconstructed image and the 
PSF is made in Fig. M. To make the structures more vis- 
ible, we used a representation similar to that of Harder 
and Chelli ( Harder and Chclli 2000| ) . The images are neg- 
atives of intensities and the representation uses threshold 
and saturation. Doing so, a strong (white) secondary max- 
imum appears in the image of P Cyg, and two lower ones 
remain far away from the core. Moreover, Fig. M clearly 
evidences that the central (white) surface of the recon- 
structed object is larger than that of the PSF (this will 
be interpreted in the next section as an effect of the enve- 
lope). It seems difficult to explain that such a structured 
image is only the result of variations of the PSF. However, 
it is also difficult to ascertain that our reconstructed im- 
age is free of any residual error (much more data would 




Fig. 7. Representation at the same scale of the recon- 
structed image of the Fig. || (left) and of the PSF (right) 
in contour plot. The white parts are due to the thresh- 
old and the saturation applied to both images: intensities 
lower than 6 % and higher than 15 % of the maximum are 
made white. North is at the top and East is at the right 
of the images. 




Fig. 8. Representations in gray levels of the residual 
blurred envelope e{r)®h{r) (top) and of the fraction "fh{r) 
of the image of 59 Cyg subtracted from P Cyg (bottom). 
The sum of these two images gives back the observed im- 
age of P Cyg. The orientation is that of Fig. n^. 



have been necessary for that). To strengthen our confi- 
dence in the fine details of the reconstructed image of the 
envelope, we have implemented a series of processing and 
tests. They are described in next sections 3.2 and 3.3, and 
in the Appendix. 
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3. 2. Deconvolution of the envelope - alone of P Cyg 

We have further processed the data assuming that P 
Cyg observed in Ha may be fairly represented by an 
unresolved star surrounded by an extended envelope. Let 
us write P Cyg as e(r) + 7<5(r), where e(r) is the envelope 
and 7 the relative intensity of the star defined by the 
Dirac function S{r). The observed image is then modeled 
as y{r) — e{r) ® h{r) + ^h{r). We have implemented a 
somewhat heuristic procedure that consists of subtracting 
"fh{r) from y{r) to obtain e(r) (g) h{r). The procedure 
has some similarity with what is made in the algorithm 
CLEAN (iHogbom 197^ , and the algorithm of Lucy 
( Lucy 1994 ). In practice, it consists of subtracting an 
appropriate fraction of the image of 59 Cyg correctly 
shifted from P Cyg. The parameters of this subtraction 
are determined to leave a smooth pattern for e(r) (g) h{r), 
with no remaining bump or hole due to the unresolved 
star. We found that a fraction 7 = 20 % of the PSF 
{h{r)), shifted of a displacement (X = - 0.3 pixels, Y = + 
1.3 pixels) was to be subtracted from P Cyg {y{r)). This 
was done using an interpolation of the images with Math- 
ematica ( Wolfram S. 1999| ). Fig. || shows the residual 
blurred envelope e(r) h{r) (top) and the fraction of the 
PSF subtracted to the image of P Cyg ^h{r) (bottom). 
This procedure, performed independently of the above 
deconvolution and parametric estimation of the star plus 
envelope, lead also to the same ratio 4 for the energy of 
the envelope relative to the central star. 



We then processed the residual blurred envelope 
e(r) ® h(r) using the same RL algorithm stopped at the 
same iteration number 245. The result is given both in 
gray levels and in contour plot with linear contour spacing 
in Fig. ^ The image of the envelope is fully consistent 
with what was obtained in the raw deconvolution of Fig. 
Q. The bright spots are all found at the same position. 
Moreover, there is a bright spot very close to the star 
clearly visible in this representation in the South- West 
direction; it was only perceptible as a small deformation 
of the central star in Fig. |[ 



3.3. Analysis of the quality of the deconvolution 

We have implemented a few more computations to check 
the quality of the results we give in this paper. The 
whole deconvolution procedure was also performed us- 
ing the Image Space Reconstruction Algorithm (ISRA) 
( Daube-Witherspoon et al. 1986| ) instead of the RL algo- 
rithm. The comparison with the image reconstructed by 
the Wiener filter leads us to the iteration 303. The recon- 
structed image is presented in Fig. |l^. The two algorithms 
gave almost the same results which reinforced us in the 
confidence we have in the deconvolution results given in 
this paper. 




Fig. 9. Image reconstruction of the envelope of P Cyg in 
gray levels (top) and in contour plot (bottom). Contours 
levels are equally spaced. North is at the top and East is 
at the right of the images. 



At that point we may conclude that the deconvolution was 
carried out taking into account the problem of noise; we 
also made use of an a priori model, by assuming that P 
Cyg was the sum of an unresolved star and an envelope. A 
significant question remains: up to what precision can we 
trust 59 Cyg as an accurate PSF for P Cyg? Several ele- 
ments enable us to give a positive answer to that question: 
the seeing conditions were similar for both observations, 
and we found the same ratio of the flux envelope/star 
(about 4) before and after deconvolution. It would have 
been very convincing to have a series of observations of P 
Cyg and a reference, eventually with different seeing con- 
ditions, and have all the results that converge towards a 
unique solution. In fact, we made an elementary test that 
consisted in dividing the long exposure of P Cyg in two 
sequences. The same operation was applied to 59 Cyg. We 
have then made a deconvolution of these four resulting im- 
ages. The results are very similar to those obtained with 
the deconvolution presented in Fig. ^ and are not repro- 
duced here. In an alternative way, we give in the appendix 
the results of a series of additional tests that tend to con- 
firm 59 Cyg as a correct PSF. 

At that point, we think that we obtained all that allowed 
the present data. Of course, futures observations, with 
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Fig. 10. Reconstructed image of P Cyg with the ISRA 
algorithm stopped at iteration 303. Contours levels are 
not equally spaced and correspond to: 100, 80, 60, 40, 
30, 15, 12.5, 9, 6.5, 5, 3.5, 2, 1.2 percent of the image's 
maximum. North is at the top and East is at the right of 
the image. 



possibly a larger telescope, will be very usefull to confirm 
our first results and to further precise the morphology of 
this object. 

4. Discussion 

The present work follows a long series of observations 
to track the history of P Cyg's giant eruptions in 1600 
and 1655. Following the first attempts in this direction 



(Feibelman 1995), it clearly appeared that the star to 



its nebula brightness ratio constitutes a challenging ob- 
stacle for the optical mapping of the nebular morphol- 
ogy. Indeed only few observations, with heterogeneous 
angular resolutions are reported in the literature. Lei- 
therer and Zickgraf ( Leitherer et al. 1987 ) first published 
the detection of P Cyg's extende d nebulosity using C CD 
imaging. Later, Johnson et al. ( Johnson et al. 1992 ) re- 
ported the detection of forbidden emission lines due to 
nitrogen enriched material at 9 arcsec. The first corona- 
graphic imaging of P Cyg from the ground was obtained 
by Barlow et al. (Barlow et al. 1994) from high resolution 
long-slit data. Barlow and co-workers discovered the pres- 
ence of a 22 arcsec nearly circular shell which presents 
bright condensations of about 2 arcsec wide, mostly in 
the North (their plate scale is of 0.255 arcsec per pixel 
where a 4.3 arcsec occulting strip was used). STScI coro- 
nagraph imaging, using a 4.8 arcsec occulting disk, shows 
almost the same features (Nota et al. 1995|). From another 



point of view P Cyg constitutes an ideal target for test- 
ing new and different imaging techniques. Among these 
there are the spectral-line image sharpening techniques 



SCASIS (|de Vos 1994|), the A MOS adaptive optics obser- 
vations ( p\4orossi et al. 1996 ), and occulting mask imag- 
ing (4 arcsec) with the new MOMI instrument for wide 
field imaging ( [O'Connor et al. 199S ). In the latter case 



the authors suggest highly asymmetrical features at 3-4pc 
from the star (7') probably associated to previous mass- 
loss events ( [Meaburn et al. 199E ). At the same time, ra- 
dio imaging now offers approximately the same panel of 
field and spatial resolution as optical imaging making the 
comparison of optical and radio maps possible. Indeed, 
sounding different scales in radio wavelength can easily be 
achieved by changing the baseline configuration of radio 
arrays (Skinner et al. 1998). Using this possibility Skin- 



ner and co-workers have compared radio maps to Bar - 
low's coronagraphic image of P Cyg ( Barlow et al. 1994 ). 
These authors claim that the emissive regions in radio and 
visible are roughly the same, although this comparison is 
further complicated by the coronagraphic mask and the 
telescope diffraction pattern. They conclude that both ra- 
dio and optical maps exhibit the same details having the 
same physical origin: i.e. dense clumps overtaken by the 
faster wind and heated by shocks. 

Low and intermediate spatial resolution images suggest 
a global spherical expanding envelope but dumpiness is 
present in each case. This trend is also present at higher 
spatial resolution of the 250km MERLIN centimetric net- 
work dSkinner et al. 1997 ) approaching the submilliarc- 
second observations of the GI2T optical interferometric 



( IVakih et al. 19971 ). 

In this context our present reconstructed image of P 
Cyg's environment in Hq presents, as expected, strongly 
clumped features within the 1 arcsec field of view (O.Ollpc 
with D=1.8Kpc). More than 6 enhanced emission clumps 
are counted with our 0.05 arcsec spatial resolution in a 
nearly 0.6 arcsec region superimposed to the unresolved 
central star. The mean size of the clumps is roughly 0.08 
arcsec which is the angular limitation of the 152cm OHP 
telescope. These results agree well both in size and mor- 
phology with MERLIN observations, at nearly the same 
resolution. The typical diameter of emitting regions for 
MERLIN is 0.4 arcsec (0.13 arcsec for the core), and 
amazingly comparable to the optical structures (Fig. g|). 
In this same figure, a North-East/South- West preferential 
axis appears in the Hq image due to the grouping of the 
clumps, the distribution being otherwise rather uniform. 
The same orientation was also pointed out in SCASIS 
observations at a lower spatial resolution ( ie Vos 1994 ). 
Note that in our reconstructed image a bright feature 
is located at 80 mas South-East of the central star. We 
can speculate on its relation to the local strong emis- 
sion discovered at 0.8 mas from the star in august 1994 
by the GI2T interferometer ( [Vakih et al. 1997| ) although 
the E-W absolute position of the latter emission was not 
given by this instrument. If this scenario holds, this po- 
sition some 3.2 years after implies a projected velocity 
around llOkm/s. Taking into account the radial veloc- 
ity and uncertainties 208±78km/s obtained by Vakili et 
al. (Vakili et al. 1997), this projected velocity is to be ex- 
pected for a clump nearly on the line of sight ejected three 
years earlier with the terminal velocity, and thus, compat- 



Chesneau O. et al.: Adaptive Optics imaging of P Cyg in Ho 



ible with the GI2T observations. Although the possible 
physical relation of 1994 blob and 1997 clump remains to 
be robustly settled, interferometric and AO imaging re- 
peated in the future, should enlight such scenario. 
At the present, only radio observations by Skinner et 
al. (Skinner et al. 1997) present confident temporal vari- 



ations. In the two 6-cm MERLIN images taken in a 40 
days interval, impressive changes were observed, corrobo- 
rated by VLA observations. In the observed region, the 
wind velocity suggests a 2 years dynamical time scale, 
which can hardly be compared to the 6-cm flux varia- 
tions. On the other hand, the recombination time scale 
for hydrogen atoms (1.210^/ne in years) is shorter, but 
not sufficient, ^160 days considering a characteristic rie 



of 2.8105c 



at 0.07 arcsec. This short time scale puts 



strong constraints on the electron density, which has to be 
four times larger than the surrounding envelope material. 
The dumpiness can explain such a time evolution if the 
structures are sufficiently small and dense, or if the shocks 
between the wind and the clumps are strong enough. The 
question is whether increasing resolution would reveal the 
same dumpiness, and if activity observed in optical and 
radio wavelengths are closely correlated. 
Some questions arise. How can these small scale clumps 
be related to the 2 arcsec ones observed in the Barlow's 
images at 3 arcsec from the star? How can they survive 
over such a long distance? Do they reappear at the loca- 
tion we detect them? 

A challenging issue is now understanding the connection 
between the different spatially resolved structures and 
their scales which needs the monitoring of the dumpi- 
ness from the star to the interstellar medium, and con- 
strains the mass ejection dynamics. Therefore, our present 
observations are the first attempt to prove that an op- 
tical monitoring of the dumpiness in inner region of P 
Cyg's mass l oss is observationa lly feasible. As previously 
pointed out ( Vakili et al. 1997 ), the intermediate regions 
of P Cyg's wind, from a few stellar radii to a few par- 
sees can be sounded by means of AO plus coronagraphic 
imaging from the ground, in relation with radio observa- 
tions. The dynamical time scale for optical interferome- 
try is around the month, but a temporal monitoring of P 
Cyg by this technique requires both higher sensitivity and 
larger numbers of baseline orientations due to the complex 
structures which occur at different scales. For larger dis- 
tances, the recombination time scale in the clumpy wind 
should produce large effects as detected in radio, and AO 
becomes the perfect technique to follow such activities. 
The brightness of P Cyg and its evolutionary time scale al- 
low the development of a multi-site and multi-wavelength 
observations campaign, using AO, optical and radio inter- 
ferometry to get a unified picture of P Cyg's environment 
physics. 




Fig. Al. Thick line:Poisson-Mandel transform P{y) of 
P{b) calculated using relation Al. Thin line: Histogram 
of the values of y{r). A good agreement between the two 
curves is observed. 



Appendix A: Tests on the PSF 

We have made a series of tests that tends to strengthen 
our confidence in the 59 Cyg's image as a correct PSF for 
P Cyg. These tests are negative tests in that sense that 
if the PSF fails to succeed one, it should be rejected as a 
bad PSF. 

Let us assume that 59 Cyg is a good unbiased representa- 
tion of the true PSF h{r). Differences between the blurred 
reconstructed image 6(r) — h{r) ® x{r) and the observed 
data y{r) must then be dominated by statistical noise fluc- 
tuations, with no bias term. In this relation x{r) is our best 
reconstructed image, as shown in Fig. ^, and h{r) the PSF 
59 Cyg of Fig. || of the body of the paper. We assume that 
the noise comes from a photo-detection process, and that 
y{r) is a realization of the Poisson process of mean b{r). 

A.l. Test 1: Poisson-Mandel transform 

The first test we have performed is a basic one, not very 
sensitive to the exact value of h{r), but that must be 
verified in any case. Let us denote y and b the values 
taken by y{r) and b{r). The total probability theorem 
(Papoulis 1984) allows us to write the unconditional prob- 
ability P{y) of y{r) as the sum of P{y/b) x P{b) for all b 
values. For the Poisson process, the conditional probabil- 
ity of y assuming b is P{y/b) — ^ exp —6 where y is an 



integer (the number of photons) and b a continuous value. 
As a consequence, P{y) is the Poisson-Mandel transform 
(iMandel 1959| [Mehta 1970[) of P(6): 



(Al) 



f°° by 

P{v) = / P(6)-rCxp-6 db 
Jo y! 



We have verified that our data correctly obeys relation [All. 
We have taken for P{b) the histogram of b{r), applied the 
above transformation and compared the result P{y) with 
the direct histogram of the values of y{r). The comparison 
is shown in Fig. Al. The results are consistent with the 
data. 
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Fig. A2. Thick line: Polar plot of the Euclidean distance 
between y{r) and bg{r) for 9 varying between and 27r. 
Thin line: Circle with a radius of the minimum of the 
Euclidean distance between y(r) and bo{r). The minimum 
of the Euclidean distance is obtained for 9 — 0. 

A. 2. Test 2: Rotations of the true PSF 

The second test we have performed was to check the cor- 
rectness of the small departures from circular symmetry 
of 59 Cyg. We use as the PSF /ie(r) an image of 59 Cyg 
rotated of an angle 9. Then the deconvolution procedure 
is carried out as previously and leads to an image X0{r). 
The blurred image hg{r) is obtained as the convolution 
of xg{r) and hg{r), and we finally compute the Euclidean 
distance between y{r) and bg{r). The results are shown in 
Fig. ^3 in a polar plot for 9 varying from to 2 tt. The 
original PSF gives the best result. If the deviation from 
circular symmetry were purely random, the goodness of 
the deconvolution would not be affected by this rotation. 
Secondary minima of the curve appear for apparent sym- 
metries of the PSF. 

Representations of the difference between y{r) and bg{r) 
are shown in Fig. A3. For the sake of clarity we have rep- 
resented the sign of the difference y{r) — hg(r). For = 0, 
we get a speckle-like pattern, roughly uniform over the 
whole image. For other values of 6, this difference shows 
large patterns that indicates regions over which hg{r) does 
not correctly matches y{r). These trends clearly show that 
hg{r) for 7^ is a biased version of h{r). 
As a conclusion, 59 Cyg passed the tests and could not be 
rejected as a bad PSF. 
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